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ABSTRACT 

We identify a physical mechanism that would have resulted in rapid, obscured growth 
of seed super-massive black- holes in galaxies at 2; > 6. Specifically, we find that the 
density at the centre of typical high redshift galaxies was at a level where the Bondi 
accretion rate implies a diffusion speed of photons that was slower than the gravi- 
tational infall velocity, resulting in photons being trapped within the accretion flow 
and advected into the black-hole. We show that there is a range of black-hole masses 
(Afbh ^ 1O'^~^M0) where the accretion flow traps radiation, corresponding to black- 
holes that were massive enough to generate a photon trapping accretion flow, but 
small enough that their Bondi radii did not exceed the isothermal scale height of self- 
gravitating gas. Under these conditions we flnd that the accretion reaches levels far 
in excess of the Eddington rate. A prediction of this scenario is that X-ray number 
counts of active galactic nuclei at z > 6 would exhibit a cutoff at the low luminosities 
corresponding to black- hole masses below ^ lO'^M©. At low redshifts we find photon 
trapping to be unimportant because it could only occur in rare low spin halos, and 
would require black-hole masses in excess of expectations from the observed black-hole 
- halo mass relation. The super-Eddington growth of ^ lO^M© seed black-holes at 
high redshift may have provided a natural acceleration towards the growth of super- 
massive black-holes at 2; '-^ 6 — 7, less than a billion years after the Big Bang. 



Key words: 

theory 



galaxies: formation, high- redshift, nuclei, quasars: general — cosmology: 



> 1 INTRODUCTION 



X 



Exceptionally bright quasars with redshifts up to 2 ~ 7 have 
been discovered (e.g. |Fan et al.||2001[ |2003[ [Mortlock et al. 
[20TT| ). These quasars are thought to be powered by the thin 



disk accretion of gas onto super-massive black-holes at the 
centres of galaxies. Their maximum (Eddington) luminos- 
ity depends on the mass of the black-hole, and the brighter 
quasars are inferred to have black-holes with masses of more 



than a few billion solar masses. Since their discovery (Fan 
et al.|200T 20031, the existence of such massive black- holes 



at z > 6 has posed a challenge to models for their forma- 
tion. This is because a ~ lO^M© black-hole accreting at the 
Eddington rate with a radiative efficiency of 10% requires 
almost the full age of the Universe at z ~ 6 to grow from 
a stellar mass seed. Many authors have therefore discussed 
solutions to this apparent mystery by including a significant 
build-up of mass through mergers ( Haiman &: Loeb||2001 |, 
collapse of low-spin systems ( Eisenstein &: Loeb||1995 1 , and 
suppression of molecular line cooling via a large Lyman- 



Werner flux (^Dijkstra ct al."2008'). Other authors have also 
attempted to explain the fast growth of black-holes at high 
redshifts based on super-massive stars ( Loeb fc Rasio|1994 



Bromm & Loeb'2003|) and more recently, quasi-stars (IVolon 
teri fc Begelman 2010l |Begelman|[20To| 



Volonteri & Rees ( 2005 1 discussed the possibility of 



super-Eddington accretion of black-hole seeds in high red- 
shift galaxies. They pointed out that seed black-holes are lo- 
cated at the centres of isothermal disks where the conditions 
for quasi-spherical Bondi accretion should be prevalent. At 
high redshift these disk centres are sufficiently dense that 
the Bondi accretion rate greatly exceeds the Eddington rate. 



Volonteri & Rees ( 2005 1 point out that this super-Eddington 



accretion provides a route by which a large fraction of the 
mass e-foldings needed to grow a super-massive black-hole 
by redshift six could be accommodated within a small frac- 
tion of the age of the Universe. However the calculations 
in Volonteri & Rees ( 2005 1 ignored feedback effects like gas 



heating, which may raise the sound-speed, and hence lower 
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the density and therefore the Bondi accretion rate. For ex- 



ample Milosavljevic et al. ( 2009 1 show that photoheating 



and radiation pressure from photoionization significantly re- 
duce the steady-state accretion rate and potentially render 
a quasi-radial accretion flow unsteady and inefficient. They 
find that the time-averaged accretion rate is always a small 
fraction of the Bondi accretion rate. Thus, the very high ac- 
cretion rates implied by the Bondi accretion in the centre 
of a high redshift isothermal disk might never be reached. 
On the other hand, if the accretion rate is sufficiently high 
that the emergent photons are trapped within the accretion 
flow, then these feedback effects cannot operate ( [Begelman] 
19791, and so the accretion rate can reach arbitrarily high 



levels. 

In this paper we flnd that at sufficiently high redshift, 
the central densities of galaxies imply Bondi accretion rates 
that exceed the rate required to trap radiation and advect it 
into a black-hole (Bcgclman 1979). Thus, we find that there 
were periods in the growth of black-holes at high redshift 
where the growth was super-Eddington and feedback mech- 
anisms could not halt the accretion flow. Our goal is not 
to make a self-consistent model for both the transport of 
material from large galactic radii and central black-hole ac- 
cretion. Rather, we note that there is a significant literature 
looking at the problem of super-Eddington accretion assum- 
ing large mass-delivery rates to the region of the black-hole, 
and identify the cosmological conditions that would provide 
sufficient accretion rates to allow this by trapping radiation. 
We begin in §[2] with a description of our simple model, be- 
fore presenting our results in § [3] We finish with a discussion 
in § |4j and conclusions in § [SJln our numerical examples. 



we adopt the standard set of cosmological parameters (Ko- 
matsu et al.pOll l, with values of f2b = 0.04, fim = 0.24 and 
Qa = 0.76 for the matter, baryon, and dark energy frac- 
tional density respectively, h = 0.73, for the dimensionless 
Hubble constant, and as = 0.82. 



2 MODEL 

The basis of this paper is a comparison between the Bondi 
accretion rate, and the accretion rate required to trap pho- 
tons within the accretion fiow. We discuss these in turn. 



Here Rd is the characteristic radius of an exponential disk 
of surface density profile 

E(r) = Eoexp(-r/i?d), (4) 

with Eo — Mdisk/27r_R|, and the disk scale length is given 
by 



-Rd = ^fvir 

V2 



(5) 



where A is the dimensionless spin parameter of the halo. In 
equation ([3|, Mdisk ~ "^dA^halo is the disk mass, rrip is the 
mass of hydrogen, and fi = 1.22 is the molecular weight 
of primordial neutral gas. At the high redshifts of interest, 
most of the virialized galactic gas is expected to cool rapidly 
and assemble into the disk. We therefore assume md = 0.17. 
The corresponding mass density is po = rripno. The virial 
radius of a halo with mass Mhaio is given by the expression 



10 



is close to unity and defined as C, = 

-31-1 



where C,(z 

[(an/!^jn)(Ac/187r2)], VL^^ = [1 + (f7A/J^m)(l + Z, , , 

Ac = 187r^-|-82d - 39d^ , and d = 0.^^-1 (see equations 22- 
25 in [Barkana fc Loeb|[200T] for more details). From equa- 
tions (|3|) and ([6|) we see that the central density no and hence 
the Bondi accretion rate MBondi scales as (1 -I- z)*, and that 
as a result accretion rates are expected to be much larger at 
high redshift. 

Dormant central black-holes are ubiquitous in local 
galaxies ( Magorrian et al. [1998 1. The masses of these super- 
massive black-holes scale with physical properties of their 
hosts (e.g. [Magorrian et al.|1998l [Merritt fc Ferrarese|2001 



Tremaine et al. 20021. However at high redshift the rela- 



tions observed in the local Universe may not be in place. 
We therefore do not impose a model for the relation be- 
tween black-hole and halo mass in this paper, and instead 
explore a range of values. Indeed, our results indicate that 
feedback, which is thought to drive the black-hole - halo- 
mass relation relation, would not be effective at early times. 
The grey curves in Figure[l]show the Bondi accretion rate as 
a function of redshift for different values of halo and black- 
hole mass. Here we assume A = 0.05 corresponding to the 
mean spin parameter for dark- matter halos ( |Mo et al.|1998 1. 



2.1 The Bondi accretion rate 



We begin with the expression for the Bondi accretion rate 
( Bondi|[T952 ) onto a central black- hole of mass Afbh 



A'/Sondi = 47rpo''Bondit'ff = 4:nno firUpr Bondi \ h^^^^ : 

V ''Bondi 

where va is the free-fall time at the Bondi radius 

_ GMbh 



(1) 



(2) 



and Ca is the sound speed, which for an isothermal gas we 
assume corresponds to a temperature of lO'^K. To evaluate 
A^Bondi we Specify the central number density of a self grav- 
itating disk ( |Schaye|2004 | 



no 



127rcii?j/imp 



(3) 



2.2 Photon trapping by the accretion flow 

If the diffusion velocity of photons at a radius r is smaller 
than the free-fall velocity of the material at radius r, then 
photons become trapped in the accretion fiow ( |Begelman| 
|1979| ). In such cases, the black-hole would be obscured. 
In this section we estimate the accretion rate required to 
achieve this photon trapping at radius r. 

The free-fall time from radius r to a smaller radius rmin 

is 



iff 



(r- 



■Uff 



(6) 



where ~ ^J~GM^^^Jr within the Bondi radius. This should 
be compared with the diffusion time of 



idiff 



_{r- 



(7) 
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Figure 1. The grey curves show the Bondi accretion rate as a function of redshift. The dark lines show the critical accretion rate for 
which the photon diffusion speed is smaller than the gravitational free-fall speed. In cases where the Bondi accretion rate is larger than 
the critical accretion rate, photons are trapped and the AGN is obscured. Three cases are considered for the vertical structure of the 
disk, which is assumed to be set by the sound speed of the gas {Left Panels), the turbulent velocity associated with a Q = 1 disk {Central 
Panels), and the virial velocity {Right Panels) respectively. In each case the Upper and Lower panels correspond to cases with halo 
masses of M = lO^M© and M = 10^°Mq. Three curves in each case correspond to black-hole masses of lO^M©, IO^'Mq and lO'^Mo. 
The cross-section was assumed to be F^ig = 100 times larger than the Thomson cross-section, the minimum radius was Fmin = 10'' times 
rg, and the spin parameter A = 0.05. 



Here the optical depth is given by 
_FsigCrT 



(8) 



the scattering cross-section of the accreting gas to the emer- 
gent radiation in units of ctt, and p is the hne-of-sight aver- 
aged density between radii Vmin and r. The density within 
the Bondi radius scales as p{r) — po{r /rBondi)~^'^ ^ yielding 



o„ ■ Duiiui r/ - mm \— u.b / ' \- 

2po— 1( ) -( j 

^min ^Bondi ^Bondi 



2p- 



-[(: 



1] for r < Tbo 



The condition t^ia > iff at radius r is satisfied for ac- 
cretion rates M > Mum where 



• p f 47rmpC 



(9) 



in which we have used the relation Meondi = 4'jTpr^VH. We 
find that the accretion rate required for photon trapping 
has a minimum value (i.e. dMum/dr = 0) at a radius of 
r — 4rmin, yielding a minimum required accretion rate for 
photon trapping of 



FsigCTT 



(10) 



Since the Eddington accretion rate at an efficiency e is 
Msdd = (47rGAfbh?Ti,p)/(ccrTe), we find 



Mlin, = 2 



mm 
Tg 



Msdd, 



(11) 



where rg = GM^oh/f? ■ 

The limiting accretion rate depends on the emergent ra- 
diation spectrum. As described below in § |3.2[ we find that 
the condition for photon trapping is more readily achieved 
for X-ray photons than for optical/UV photons. We there- 
fore frame our discussion around trapping of optical/UV 
photons. In the case of optical/UV photons the cross-section 
would be dominated by dust at radii larger than a sublima- 
tion radius rmin = ^sub, beyond which the opacity of the 
gas to optical/UV photons could be larger than the Thom- 
son opacity by as much as two or three orders of magnitude 
due to the presence of dust (Laor & Draine, 1993). The ex- 



istence of dust can be questioned for high-redshift galaxies 
(e.g. Bouwens et al. 20101, however high metallicities are 



inferred from the broad emission lines of all quasars out to 
z = 7.1 and so metal enrichment (due to star formation) is 
known to precede the growth of the black-hole in the galac- 
tic nuclei of interest ( Hamann fc Simon|2010 1 . On the other 
hand, the diffusion time may be lessened if the gas is clumpy, 
corresponding to a lower effective opacity. We therefore take 
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a typical value for the cross-section that is F^ig = 100 times 
larger than the Thomson cross-section (additional cases are 
presented below in ; 



a typical value of Fa, 



2.3 1. Defining rmin = F^inTg we assume 
in = 10"^ to describe the sublimation ra- 



dius ( |Netzer fc Laor|1993 |. 

Limiting accretion rates corresponding to these default 
values are plotted in the left hand panels of Figure [l] (dark 
lines) for the case where the vertical structure of the disk 
is set by the sound speed of the gas. The Upper and Lower 
panels correspond to cases of halo masses of A/haio = W^Mq 
and Mhaio ~ 1O^''M0. The larger value corresponds to the 
lower end of the inferred halo masses of the Lyman-break 
population at 2 > 6 ( [Trenti et al.|2010| . 



3 PHOTON TRAPPING 

In cases where the Bondi accretion rate is larger than the 
critical accretion rate, rest frame optical/UV photons are 
trapped and the active galactic nucleus (AGN) can be ob- 
scured. Thus the cross-over of the limiting and Bondi ac- 
cretion rate curves in Figure [T] represents the redshift be- 
yond which accretion traps radiation for an AGN powered 
by the particular black-hole masses shown. If the vertical 
disk structure is set by the sound speed, we find that a 
lO^M© black-hole in a 1O^°M0 halo will result in trapped 
radiation at 2; > 4. In this section we discuss the range of 
halo and black-hole masses that result in photon trapping. 



2.3 Disk structure 

Before proceeding we first note that we have so far as- 
sumed the thickness of the gaseous disk at the centre of 
the galaxy to be set by the sound speed of gas at lO^K. 



However from equation ( 13 1 we see a strong dependence on 



the value for the effective sound speed. For example, in tur- 
bulent disks, the turbulent velocity replaces the isothermal 
velocity in determining the effective sound speed. Recently, 



Genzel et al. (20101 inferred a Toomre's Q value of Q = 1 



in ULIRGS, implying a value for the turbulent velocity of 
GT,/n = ^GSor/vr, where = GT,onr^/r, yielding 



ct 



Q — v/r = ^ GY^Q-K jr. Evaluating at the Bondi radius, we 
get 



Ct 



(12) 



This value of ct is the maximum value possible for a disk 
at large radius (as a higher ct would imply an unphysical 
disk with h > r). Therefore, for a Q = 1 disk the turbu- 
lent velocity should decrease towards small r. At sufficiently 
small radii, ct < Cs, implying that the minimum thickness 
of the disk at its centre is set by the sound speed in the 
gas. Thus, although the turbulent velocity can be signifi- 
cant at large radii in the disk, we find that when evaluated 
at the Bondi radius, the turbulent velocity is smaller than 
the sound speed unless the halo is very massive and at high 
redshift (> lO" at z > 4). 

There may also be sources of energy injection that heat 
the gas to temperatures that correspond to velocities much 
larger than ct (or Cs), although in order to remain bound the 
velocity of the gas must be smaller than the virial velocity 
of the halo. We therefore show two cases in addition to Cs ~ 
lOkm/s, namely ct and Cv ~ fvvii where / ~ 0.5 in order 
to bracket the range of interest, with corresponding limiting 
accretion rates plotted in the central and right hand panels 
of Figure [1] (dark lines). The results are almost unchanged 
where ct is used to set the disk height. This is because we 
find that for the halo masses, black- hole masses and redshifts 
of interest, the sound speed is smaller than the turbulent 
velocity when evaluated at the Bondi-radius. However if the 
disk height is set by the maximum velocity Vc, then larger 
black-holes and smaller halos are needed in order for the 
Bondi accretion rate to exceed the limiting rate, and at 2 ~ 6 
black-hole accretion would not be obscured in this case. In 
the remainder of this paper we restrict our attention to the 
case of an isothermal disk. 



3.1 Scaling relations 

The conditions for the halo mass, black-hole mass, and red- 
shift that conspire to provide accretion rates that trap the 
optical/UV radiation can be obtained by combining equa- 
tions |T| and (10 1. Evaluating we find 

2/3 



A/bondi 



37.5 



F ■ 

^ mm 



( Mbh 
VIOSMq 
1 



f A/halo 
VIOIOM0 



lOkm/s 



100 / 



1 + 2 
7 

A 
005 



0.17 



(13) 



This expression makes explicit the dependencies that lead 
to photons being more easily trapped (i.e. Mbondi/Mum > 
1), namely larger black- holes in larger halos, and at higher 
redshift. We note that the condition for photon trapping has 
a much steeper dependence on redshift than obscuration via 
dust absorption, which scales as (1-1-2)^. 

This critical redshift at which the Bondi accretion rate 
exceeds the critical rate for photon trapping is shown graph- 
ically in Figure [2] Here we plot contours of the redshift 
at which the Bondi accretion rate becomes larger than the 
critical accretion rate, as a function of halo and black-hole 
mass. This figure shows the mass combinations that give 
a cross over from AGN in which radiation can escape to 
those in which the radiation is trapped. Based on equa- 



tion ( 13 1 we plot contours for the parameter combinations 

(Fsig/l#)(i^mm/10^)-' = 1, (^'sig/10')(F„,in/103)-l = 0.3 

and (Fsig/10^)(Fmin/10^)"^ = 3 in the Left, Central and 
Right panels of Figure[2] Smaller (larger) values of this com- 
bination lead to larger (smaller) inflow rates being needed 
to trap the radiation, and so larger (smaller) black-holes are 
needed to trap the radiation at flxed halo mass and redshift 
(see equation fTsl . 



From equation (111, the relation between the limiting 



rate for photon trapping and the Eddington rate can be 
illustrated via the condition 



Mlim 
A^Edd 



102 



Fm 



103 



(14) 



An accretion rate that leads to photon trapping is likely to 
be in excess of the Eddington rate, and hence associated 
with a very rapid build-up of black- hole mass. 



3.2 Trapping of X-rays 

Up until now we have found that the critical rate at which 
the accretion flow traps optical/UV photons can be reached 
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Figure 2. Contours of the redshift at which the Bondi accretion rate becomes larger than the critical accretion rate as a function 
of black-hole and halo mass. We show contours for cases where the vertical structure of the disk is assumed to be set by the sound 
speed of the gas. In the Left, Central and Right panels we show contours for parameter combinations (-Fsig/10^)(-Fmin/10^)~^ = 1, 
(Fsig/10^)(i^inin/10'^)~'^ = 0.3 and {i^sig/10^)(-Fjni,i/10'')~^ = 3. For illustration, the grey regions show the portion of parameter space 
at z ~ 6 that does not result in photon trapping accretion flows. The example assumes A = 0.05 and = 0.17. 



in the centres of high redshift galaxies owing to gas with 
large opacity (i.e. Fsig ^ 1) beyond a radius Fmin»'g. The re- 



sults of Laor & Draine ( 1993 1 indicate that while Faig ^ 1 is 



expected for UV photons, the X-ray component of the spec- 
trum will see an opacity to the inflowing gas that is set by 
the Thomson cross-section (i.e. Fsig^x = !)■ Moreover, if the 
X-ray photon component of the spectrum were not trapped, 
it could halt the accretion flow via radiation pressure if it 
exceeded the Eddington rate by itself, thus preventing trap- 
ping of the optical/UV radiation. On the other hand, X-rays 
see this Thomson opacity at radii much smaller than the 
sublimation radius, with Fmin,x ~ 10. 

To illustrate the importance of X-rays in this context 



we modify equation ( 14 1 describing the relation between 



the limiting rate for optical/UV photon trapping and the 
Eddington rate that is specific to the X-ray portion of the 
spectrum (A^Edd.x). We define Fx to be the fraction of the 
luminosity in X-rays, which based on the spectrum in |Elvis| 
|et al.| fTgtit takes a value of Fx ~ 10%. We therefore find 



0.2 



Bdd.X 



0.1 



102 



10^ 



-Fx 
0.1 



(15) 



Here Fgig ^ 1 and Fnin ^ 1 correspond to the val- 



ues for UV photon trapping. Equation ( 15 1 shows that in 



order for X-rays to exceed the Eddington limit, the pa- 
rameters corresponding to optical/UV trapping need to be 
(Fsig/10^)(Fmin/10^)"^ > 5. At these large accretion rates 
we find that when optical/UV photons are trapped, so too 
are the X-rays unless the minimum radius at which the X- 
rays encounter opacity (rmin.x = Fmiii,x''g) is described by 
a value of Fmin.x > 50. To see this we note that the ratio of 
the Bondi accretion rate to the rate needed to trap X-rays 
is 



Mb 



-Fmin.X 



1 



(16) 



Since we expect -Fmin.x ~ 6, corresponding to the innermost 
stable circular orbit, we expect that an accretion flow which 
traps the optical/UV radiation will also trap the X-rays. 



As a result we do not consider the effect of X-rays on the 
accretion flow for the remainer of this paper. 



3.3 Critical spin parameter 



Equation ( 13 1 shows that the trapping of photons is very 



sensitive to the value of the spin parameter A which governs 
the density of the galactic disk. In particular, photons may 
be more easily trapped within disks of low spin parameter. 
The assembly of dark matter halos leads to a distribution of 
spin parameters. To illustrate the parameter space we recast 



equation ( 13 1. Photon trapping at all wavelengths occurs for 
spin parameters A < Ac where 



Ac ~ 0.12 



Mb 



1O5M0 



1/4 



/ Mhalo 



VIOIOM0 



1/6 



l + Z 



103 



-1/4 



-5/4 



F ■ 

^ sig 
102 



rrid \i/2 
017/ 

1/4 

, (17) 



lOkm/s^ 

In the left panel of Figure |3] we plot contours of Ac as a 
function of black-hole and halo mass at z — 6. The example 
assumes (Fsig/10^)(Fmin/10^)"^ = 1, and md = 0.17. The 
mean spin parameter A = 0.05 is shown in black, and the 
grey shaded region illustrates the region of parameter space 
where trapping is not possible for a mean disk. Figure [S] 
shows that emergent photons from black-holes with Afbh 
1O^~'*M0 will be trapped at all wavelengths in the centres 
of high redshift galaxies (~ 1O*~^''M0) that are at the mean 
of the spin-parameter distribution . 



3.4 Breakout of the Bondi radius 

In the previous sections we have illustrated that the large 
densities expected at the centres of high redshift galaxies 
lead to conditions where accretion rates may be sufficient 
to trap photons within the accretion flow. These calcula- 
tions are based on the density at the centre of a pressure 
supported self-gravitating disk. In this section, we point out 
that the calculation applies only to black-hole masses for 
which the Bondi radius is smaller than the scale-height at 
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log(My,a\o/MQ) log(My^a^o/MQ) log(My^a^^/MQ) 

Figure 3. Regions of black-hole - halo mass parameter space which result in photon trapping. Left: Contours of Ac as a function of 
black-hole and halo mass. The grey shading shows the region of parameter space which does not result in photon trapping for the mean 
disk, because the critical spin parameter is smaller than A = 0.05. Center: Contours of Ac,b as a function of black-hole and halo mass. 
The grey shading shows the region of parameter space which does not result in photon trapping for the mean disk, because the Bondi 
radius is larger than the scale height at the disk centre. Right: Contours of the probability that accretion will result in photon trapping, 
corresponding to the spin parameter A lying between the two critical values, i.e. A^ b < < -^c To calculate the probability we assume 
the distribution is Gaussian in the natural logarithm In A, with variance a\ = 0.5 and a mean at A = 0.05. Contours arc shown for 
P = 1%, P = 10% and P = 50%. The examples assume z = 6, (-Faig/10^)(-Fmin/108)-i = 1, and = 0.17. 



the disk centre. Moreover, we note that the trapping of ra- 
diation cannot be reahsed once the black-hole grows suffi- 
ciently that its Bondi radius exceeds the scale height. 

We again utilise the self-gravitating disk model. The 
scale height at the disk centre is 



20 



ttGEo' 



(18) 



This expression ignores the gravitational contribution from 
the black-hole which would serve to reduce the scale height, 
and is valid when rBondi < 20 • Photon trapping and obscura- 
tion are only possible in this regime. We therefore calculate 
the ratio of Bondi radius to the central disk scale height as 



fBondi 
20 



(19) 



the critical value of spin parameter Ac,b at which rsondi ~ 20 



0.08 



A4 



lOkm/s 



10iOM( 
1/2 



(21) 



In the central panel of Figure |3] we plot contours of Ac,b 
as a function of black-hole and halo mass at 2 = 6. The 
example again assumes (-Fsig/10^)(-Fmin/10'^)^^ — 1, and 
nid ~ 0.17. The mean spin parameter A = 0.05 is shown 
in black, and the grey shading illustrates the region of pa- 
rameter space where photon trapping is not possible for the 
mean disk because the Bondi radius is larger than the scale- 
height. Figurejsjshows that black-holes with Mbh ^ 1O'*'^M0 
within high redshift galaxies (~ lO^'^M©) having the mean 
of the spin-parameter cannot form a photon trapping accre- 
tion flow. 



Putting in characteristic values we obtain 



''Bondi 

20 



= 2. 



VlOSMoj VIOIOM0 



1/3 



V0.17/ 



Cs 

lOkm/s 



1 + 2^ 

7 

4 -2 

0.05 



.(20) 



Thus, within M ~ IO^'^Mq halos at 2 ~ 6, black-holes 
in excess of Mbh ~ lO^M© have Bondi radii which are larger 
than the disk scale-height and so may not be obscured. We 
note that the ratio is very sensitive to the value of the sound 
speed, with a value larger than the fiducial lOkms"^ sig- 
nificantly reducing the ratio, allowing larger black-holes to 
accrete in the photon trapping mode. A small gas fraction 
also reduces the ratio. However the breakout of the Bondi 
radius implies that high mass black-holes such as those ob- 
served within the SDSS quasars (e.g. |Fan et al.|2001| |2003| ) 
could not have their emergent radiation trapped. 

To better understand the constraints imposed on black- 
hole masses where photon trapping can occur, we evaluate 



3.5 When and where could photon trapping 
occur? 

The right panel of Figure [3] shows the probabilities that the 
spin parameter A lies between the two critical values needed 
for the conditions of i) photon trapping, and ii) a Bondi 
radius that is contained within the disk scale height (i.e. 
Ac,B < A < Ac). To calculate this probability 



P = 



27r(TA 



exp 



(InA-lnA)^ 



2al 



dX 



(22) 



we assume the distribution of spin parameters is Gaussian 
in the natural log arithm In A, wit h variance a\ = 0.5 and a 
mean at A = 0.05 ( |Mo et al.|l998[ ). The probability is P = 
if Ac < Ac,B, indicating a black-hole - halo mass combina- 
tion that cannot produce a photon trapping accretion flow. 
Contours are shown that represent black-hole - halo mass 
combinations for which P = 1%, P = 10% and P = 50% of 
disks would have densities that result in photon trapping. As 
before, the example assumes (-Fsig/10^)(-Fmin/10^)~^ — 1, 
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and TTid = 0.17. The mean disk at 2 ~ 6 has a central 
density that leads to photon trapping for black-hole masses 
up to Mbh W^Mq within halos up to Mhaio ~ IO^Mq. 
However ~ 10% of disks are dense enough that photon trap- 
ping will occur for black- hole masses up to Alhh ~ 10^ Mq 
within larger halos up to Mhaio ~ IO^^Mq. Thus, we would 
expect photon trapping to be common in galaxies hosting 
Mhh ^ IO^Mq black-holes at z ~ 6. Moreover, since from 



equation ( 13 1 we see that this growth is super-Eddington, we 
find that photon trapping provides a mechanism by which 
rapid black-hole growth could proceed at high redshift, help- 
ing to explain how super-massive black-holes grow less than 
a billion years after the Big-Bang. 

In Figure |4] we explore how the conclusions regarding 
the black-hole and halo mass-ranges that produce photon 
trapping accretion flows are effected by redshift and the 
parameter combination (-Fsig/10^)(-Fmin/10'')~^. We choose 
three values of redshift, z = 1, z = 6 and z — 10. For 
z = 6 and z = 10 we choose rud = 0.17. However at later 
times we expect that gas is less plentiful and so assume 
nid = 0.025 at z ~ 1, corresponding to a comparison with 



low redshift disks (Mo et al. 19981. These cases are shown 



in the upper central and lower rows respectively. In each 
case we show examples with (Kjig/10^)(Fmin/10^)~^ = 0.3, 
(F.ig/10')(F„,in/10^)-^ = 1 and (Fsig/lO'jCF^in/lO^*)-! = 
3 {Left, Central and Right panels respectively). In all panels 
we show contours for the probabilities that the spin parame- 
ter A lies between the two critical values needed for the con- 
ditions of i) photon trapping, and ii) a Bondi radius that is 
contained within the disk scale height (i.e. Ac,b < A < Ac). 
As before contours are shown that represent black-hole - 
halo mass combinations for which P — 1%, P = 10% and 
P = 50% of disks would have central densities that lead to 
photon trapping. 

We flnd that larger values of (Fsig/10^)(Fmin/10^)*^ 
lead to more massive black-holes in smaller halos having 
their emergent radiation trapped. At redshift z ~ 6 — 10 
we find that black- holes with masses up to Mbh ~ 10^ M© 
have their radiation trapped in 10%-50% of cases. At lower 
redshift, black-holes with masses of Mbh ~ IO^ '^Mq could 
have radiation trapped in up to 10% of cases, but only if 
the halo mass and (Fsig/10^)(Fmin/10^)"^ are large. Thus, 
photon trapping is likely to be a phenomenon dominated 
by low mass (Mbh J; IQ^VIq) black-holes in high redshift 
{z > 6) galaxies (Mhaio ~ 10'°Mq). 



3.6 Comparison with Mbh — Mhaio models 

As noted in the introduction, the relations observed between 
black-hole mass and galaxy properties in the local Universe 
may not be in place at high redshift. For this reason we have 
not imposed a model for the relation between black-hole and 
halo mass in this paper, and instead have explored a range 
of values. However it is interesting to compare the range of 
black-hole masses found to be accreting in the photon trap- 
ping mode with expectations of the simple models relating 
black-hole and halo masses that have been successful in de- 
scribing some of the properties of high redshift quasars. 



Motivated by local observations ( jFerrarese 20021, we 
consider a model in which the central black-hole mass is 
correlated with the halo circular velocity. This scenario is 



quasars out to z; ~ 3 and demonstrated that the relation be- 
tween black-hole masses and the stellar velocity dispersion 
does not evolve with redshift. This is expected if the mass 
of the black-hole is determined by the depth of the gravi- 
tational potential well in which it resides, as would be the 
case if growth is regulated by feedback from quasar outflows 



(e.g. |Silk fc Rees]p98| [Wyithe fc Loeb||2()03| ). Expressing 
the halo virial velocity, Vc, in terms of the halo mass, Mhaio, 
and redshift, z, the redshift dependent relation between the 
super-massive black-hole and halo masses may be written as 



Mbh = EbhMh 



f Mhaio 

VIQIOMq 



2/3 



[C(^)] 



5/6 



l + Z 
7 



5/2 



(23) 

The normalising constant in this relation has an observeqj 



value of ebh ~ 10 based on calibration of equation (23 1 
(Ferrarese 20021, where we take the underlying 



at z 



assumption that the halo mass proflle resembles a Singular 
Isothermal Sphere. In |Wyithe fc Loeb| ( |2005[ ) this model was 
shown to be consistent with the clustering and luminosity 



function data of the 2dF quasar redshift survey ( Croom et al 
I2OO2I1. 



Interestingly, the predicted value for super-massive 
black-hole masses in galaxy mass halos (as shown by the 
thick grey curve in Figure |4| is comparable to the range 
at which we would expect photon trapping in galaxies with 
z > 6. Thus, prior to obtaining masses where the self regulat- 
ing mechanisms thought to be responsible for the relations 
between black-hole mass and halo properties take effect we 
would expect black-holes in high redshift galaxies to reach 
the level where super-Eddington accretion would become a 
natural part of their evolution. However at low redshift, Fig- 
ure|4]shows that photon trapping requires black-hole masses 
that are greater than the observed black-hole - halo mass 
relation. Thus, we do not expect expect photon trapping 
within low redshift galaxies. 



4 DISCUSSION 

The flndings in this paper point to some potentially im- 
portant implications for the growth of high redshift super- 
massive black- holes. Recent simulations by [lI] ( |2011| | show 
that self gravity overcomes radiative feedback and that 
accretion onto intermediate mass black-holes reaches the 
Eddington rate. However, at high enough accretion rates, 
a spherical inflow is not subject to the Eddington limit 
(though the emergent radiation is), and the Bondi accre- 
tion rate can take arbitrarily large values in high density 
galactic centres ( Begelman|1979 1. The conditions for spher- 



ical accretion may only rarely be realised at the centres of 
high redshift galaxies. On the other hand, there is a class of 



slim disk models (e.g. 


Abramowicz et al.| 


1988, Chen et al. 


1995 


lOhsuga et al.|2005||Watarai|2006||Ohsuga & Mineshige 


2011 


1, for which super-Eddington accretion flows are possi- 



ble. These models include a large viscosity, and so rapid 
transport of material through the disk within a small multi- 
ple of the free fall time (e.g. | Wat"Mai|2006| ) , and are optically 



supported by the results of Shields et al. ( 2003 1 who studied 



^ We have used the normalization derived by | |Ferrarese||2002[ | 
under the simplifying assumption that the virial velocity of the 
halo represents its circular velocity. 
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Figure 4. Contours of the probability in black-hole - halo mass parameter space that accretion will result in photon trapping, 
corresponding to the spin parameter A lying between the two critical values, i.e. Ac,b < A < Ac. Contours are shown for P = 1%, 
P = 10% and P = 50%. Examples are shown for three values of redshift, z = 1, z = 6 and z = 10 {Top to Bottom). For z = 6 and 
z = 10 we choose = 0.17. We assume m(j = 0.025 at z ~ 1. In each case we show examples with {i^sig/10'^)(-Fmin/10^)~"'^ = 0.3, 
(Fsig/102)(Fmin/10^)^i = 1 and (Fsig/102)(Fmin/103)"^ = 3 (from Left to Right). The thick grey curve is the predicted value for 
black-hole mass halo mass relation (equation |23|l. 



thick and advection dominated ( Chen et al.||1995 K Photon 
trapping plays an important role within these disks, even 
those with complex flows, leading to very inefficient energy 
conversion and a luminosity that is independent of accretion 
rate (e.g. Ohsuga et al.|[2005 K Simulations indicate that at 
sufficient densities, the mass accretion rate can reach hun- 
dreds of Eddington, although the photon luminosity does 
not (e.g. [Ohsuga &: Mineshige|2011[ ). These super-Eddington 
accretion disks will generate outflows (which cannot be gen- 
erated in spherical accretion, Begelman|1979 1 that are likely 
to be channelled along the poles, and so will not suppress 
the accretion along the equatorial plane. The calculations 
presented in this paper provide the boundary condition for 



super-Eddington flows. These rapid accretion events may 
provide the seeds for super-massive black-hole growth. 

A plausible scenario therefore includes two episodes of 
black-hole growth. Initially the accretion may have been 
spherical or via a slim disk with high viscosity, so that the 
Eddington rate was greatly exceeded by the accretion rate. 
We find that radiation would have been trapped and ad- 
vected into the central black-hole by the very large Bondi ac- 
cretion rates at high redshift. This phase of growth would be 
obscured at both optical/UV and X-ray wavelengths. Once 
the Bondi radius became larger than the scale-height of the 
galaxy disk, the accretion rate dropped, and higher angular 
momentum gas would have settled into a thin disk accretion 
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mode, in which the accretion time was much longer than the 
free-fall time, allowing the radiation to escape. For a pressure 
supported disk, we have found that the Bondi radius is ex- 
pected to exceed the scale-height of the disk for black-holes 
in excess of ~ 10^ Mq. The luminous quasars discovered at 
z >6 with black-hole masses in excess of IO^Mq are there- 
fore thought to be shining in this later mode, un-obscured 
in the optical and with accretion rates close to, but smaller 
than Eddington. However a portion of their prior growth, 
when the black-hole mass was ~ 1O*~^M0 would have been 
in the photon trapping mode. 

Super-Eddington accretion rates arising from the large 
densities in the cores of high redshift galaxies were previ- 



in a fraction of galaxies rather than in all galaxies as implied 



ously considered by Volonteri & Rees (20051. As noted in 
the introduction, the calculations in | Volonteri fc Rees| ( |2005l ) 
neglected feedback effects like gas heating, which may lower 
the the Bondi accretion rate (e.g. Milosavljevic et al.|2009 |. 
In this paper we note that if the accretion rate is sufficiently 
high that the emergent photons are trapped within the ac- 
cretion flow, then feedback effects cannot operate (Begel- 
| man|[T979| . Our model also makes two different assump- 
tions which modify our conclusions relative to [Volonteri fc| 
Reesl pOOSl. Firstly, IVolonteri & Reesl (120051) assume that 



once the gas is enriched, metal line cooling allows the gas 
to cool to temperatures much lower than 10*K, so that it 
fragments to form stars and the super-Eddington accretion 
episode is ended. This assumption was necessary in order 
that super-Eddington accretion not lead to black-hole den- 
sities in excess of those observed. However, we expect that 
fragmentation and associated star-formation will reheat the 
gas via radiative and supernova feedback so that the pres- 
sure support at the centre of the disk is maintained and 
accretion can continue. This scenario is supported by the 



observations of Genzel et al. (20101 which imply significant 
turbulent velocities in high redshift galaxies. Rather than 
make an arbitrary assumption that growth stops when the 
size of the accretion disk grows by a factor of five, we instead 
assume that the super-Eddington photon trapping accretion 
mode would be regulated by the time when the Bondi radius 
exceeds the scale height of the pressure supported disk. We 
find that this condition prevents super-Eddington accretion 
onto high mass black-holes. 

4.1 Obscured accretion 

Our results have some relevance to the recent discussion 
surrounding obscured accretion in high redshift galaxies. 
While luminous optical quasars in the most massive ha- 
los (> W^^Mq) dominate the observations of high redshift 



super-massive black-holes, Treister et al. ( 2011 1 recently pre- 
sented evidence that most of the black-hole accretion at 
> 6 is actually optically obscured, and in galaxies below 
halo masses of ~ 1O^°""M0. Since high redshift Lyman- 
break galaxies are thought not to be dusty (e.g. 



Bouwens 



et al. 20101, our results might have provided a mechanism 



by which the AGN can be obscured even in the absence of 
a large dusty component. However our results do not sup- 
port a photon trapping explanation for this result. Firstly, 
X-ray photons would be more strongly trapped by the ac- 
cretion flow than the UV photons, indicating that we would 
not expect to observe X-rays without optical detection. In 
addition, we also find that photon trapping is only expected 



by Treister et al. (20111 



While other mechanisms may lead to buried accretion, 
the findings of Treister et al. (2011 1 have also been disputed 
by a number of authors ( Fiore et al.|201l" Cowie et al.|2011 



Willott||2011 1, who do not observe the same level of X-ray 



emission. These authors find limits on the average X-ray 
luminosity in the rest frame 0.5-2keV band of Lo.5-2 < 4 x 
lO*'^ erg/s for z ~ 6.5 dropout galaxies. This luminosity can 
be related to black-hole mass as 



3 X 10* 



Mb 



IO^Mg 



r;ergs 



(24) 



where r; = 1 is the fraction of the Eddington accretion rate, 
and we have assumed the spectral energy distribution of 



Elvis et al. (19941. Observed luminosities must have 77 < 1 



even if the accretion rate is super-Eddington (Begelman 



19791. The observed luminosity limit therefore corresponds 
to observed black-hole masses of A'/bh < 1.3 X IO^Mq. We do 
not find that black- holes with masses Afbh -Sj 10"" M© produce 
photon trapping accretion flows, and so photon trapping 
does not explain the lack of observed X-ray sources among 
the z ~ 6.5 d ropouts ( |Fiore et aL]|2011| [Cowie et ar||2011 



Willott|2011 1 . Since the stacked observations in these stud- 
ies are based on only ~ 10^ galaxies, this lack of detection 
could follow from the low-duty cycle of AGN (which is likely 
a few percent, Wyithe & Loeb 20031. However our results 



do suggest that 90%-100% of disks with black-holes below 
this mass would be in the photon trapping mode. Thus, we 
would expect that deeper and wider held X-ray observations 
using future X-ray observatories, to display a cut-off in the 
X-ray luminosity function at about Lo.5_2 ~ 3x 10'*^ergs~^. 
Conversely, the discovery of X-ray AGN with luminosities an 
order of magnitude lower than current limits would there- 
fore rule out photon trapping accretion as a mechanism for 
rapid growth of early black-holes. Finally, we note that since 
super-Eddington accretion at high redshift is obscured at 
both optical and X-ray wavelengths, rapid growth of seed 
black-holes could not provide a significant source of X-ray 
for reionization of the IGM ( Volonteri fc Rees|2005 l. 



4.2 Seed black-hole growth 

The photon trapping mode is likely to be important for the 
rapid growth of seed super-massive black-holes with masses 
of ~ 1O''~^M0. Because the Bondi accretion rates in these 
high redshift galaxies could be orders of magnitude larger 
than the Eddington rate, the photon trapping mechanism 
helps alleviate the difficulty of growing super-massive black- 
holes of more than a billion solar masses (corresponding to 



the most distant known quasars, Mortlock et al. 2011 1 within 



the first billion years of the Universe's age. This point was 



made in detail in Volonteri & Rees (20051. To illustrate we 



note that accretion at the Eddington rate (with e — 0.1) 
leads to an e-folding time of t = 4 x lO'^ years. Assuming that 
a black-hole accretes with a duty cycle of unity, the number 
of e-folding times available by 2; ~ 7 is therefore ~ 20. This 
should be compared with the 20 e-folds needed to grow a 
IM0 black-hole seed up to a mass of IO^Mq. Thus there 
is only just enough time during the age of the Universe at 
z > 6 for a stellar black-hole seed to grow to a super-massive 
black-hole. We suggest that the period of obscured growth 
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in some galaxies would provide a path toward growing these 
super-massive black-holes. 



5 CONCLUSION 

In this paper we have determined the cosmological regime 
in which photons produced through accretion onto a cen- 
tral black-hole are trapped by infalling material, so that i) 
radiation feedback on the infall of gas outside the dust sub- 
limation radius is suppressed, allowing accretion rates far 
in excess of the Eddington limit, u) AGN appear obscured, 
and in) the growth time of black-hole is short. Specifically 
we find that a large fraction of galaxies at 2: > 6 with 
masses up to those of the observed Ly-break population 
(halo masses of ~ 1O^~^^M0) exhibit Bondi-accretion rates 
onto Mbh ~ 1O^~^M0 black-holes that are sufficiently high 
to trap the resulting rest frame optical/UV/X-ray radiation. 
The obscuration due to photon trapping is found only to 
occur for black-holes with masses up to ~ ICM© because 
larger black-holes have Bondi radii that exceed the scale 
height of the disk from which they accrete gas, so that a 
photon trapping mechanism cannot operate. As a result, we 
find a natural distinction between obscured, photon trapping 
accretion onto ~ lO^M© black-holes in galaxies of halo mass 
< IO^^Mq , and the luminous accretion seen in the brightest 
quasars with black-hole masses of ~ 1O®~^M0 within halos 
of mass ~ 1O'^^~'^^M0. At lower redshift photon trapping 
requires black-holes that are larger than expected from the 
black-hole - halo mass relation, and so is not expected to be 
observed. Our results indicate that super-Eddington accre- 
tion of mass to form seed black-holes of ~ IO^Mq provided 
a mechanism by which super-massive black-holes were able 
to form prior to ^; ~ 6. 
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